Electromagnetic propulsion and separation by chirality of nanoparticles in liquids 
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We introduce a new mechanism for the propulsion and separation by chirality of small ferromag- 
netic particles suspended in a liquid. Under the action of a uniform d.c. magnetic field H and an 
a.c. electric field E isomers with opposite chirality move in opposite directions. Such a mechanism 
could have a significant impact on a wide range of emerging technologies [J]. The component of the 
chiral velocity that is odd in H is found to be proportional to the intrinsic orbital and spin angular 
momentum of the magnetized electrons. This effect arises because a ferromagnetic particle responds 
to the applied torque as a small gyroscope. 

PACS numbers: 47.63. mf, 87.80. Ek, 87.50.ch 
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Recent years have witnessed an explosion of interest 
in the fabrication of nanoscale objects 041] an d their 
propulsion in liquids. Numerous mechanisms have been 
proposed to achieve this propulsion ranging from elec- 
trophoresis for platinum rods to beating of flexible 
magnetic rods resembling flagella 0. Some mechanisms 
take advantage of the lack of center-symmetry of the par- 
ticles and propel them under the action of rotating ex- 
ternal fields 

In this article we introduce a new mechanism for the 
propulsion and separation by chirality of small ferromag- 
netic isomers, which are suspended in a liquid and en- 
dowed with a frozen-in magnetic moment M and electric 
dipole moment d. The separation by chirality is induced 
by applying a uniform linearly polarized a.c. electric field 
E and a d.c. magnetic field H. We note that in a given 
particle the direction of the electric dipole moment d is 
unique. In contrast, the magnetic moment M can acquire 
different orientations with respect to the crystalline axes. 
Here we consider particles with an easy axis magnetic 
anisotropy in which states with the two possible mag- 
netization values, M and — M, are realized with equal 
probabilities. We consider the general case where d and 
M are non collinear, and so are E and H. In what follows, 
we will be interested in the chiral velocity V c ^ = Vq — Vq 
of the particles averaged over different directions of the 
magnetization. Here Vq and Vq are the velocities of the 
right and left-handed particles, respectively. 

Since H is an axial vector and E a polar vector the 
time-averaged chiral velocity of the particle must have 
the following form 

V ch = <7ih + er 2 (e • h) e + er 3 (e • h)[e x h]. (1) 

Here h and e are unit vectors in the direction of the 
magnetic and electric fields respectively. The coefficients 
<7j have opposite signs for isomers of opposite chirality. 
Thus, particles of opposite chirality will move in opposite 
directions. 



It is important to note that a ferromagnetic particle 
possesses an intrinsic angular momentum associated with 
the magnetization, L e = M/7 [lOj . The gyromagnetic 
ratio 7 can be estimated as 7 ~ e/mc, where m and e 
are the electron mass and charge, and c is the speed of 
light. The value of L e is relatively small, and it is usually 
neglected in studies of the dynamics of small ferromag- 
netic particles. However the existence of L e ^ means 
that the particle responds to external torques as a small 
gyroscope. This leads to new physical effects. 

We show below that the first two terms in Eq. (JlJ arc 



proportional to L e 



(the third term, on the other 



hand, remains finite even as L e — > 0). This can be seen 
from the following consideration. In the approximation 
where L e = the magnetic moment affects the particle 
motion only via the external torque 



r = MxH + dxE. 



(2) 



Thus, in this approximation all terms in the velocity of 
a particle that are odd in H should also be odd with 
respect to M. Consequently, these terms vanish upon 
averaging over different realizations of M. On the other 
hand, since L e oc M the chiral current can contain terms 
which arc odd in M and linear in L e , which do not vanish 
upon averaging over different realizations of the particle 
magnetization M. 

The motion of small particles in a dilute suspension 
can be described using the formalism of low Reynolds 
number hydrodynamics (llT | in which the external forces, 
F, and torques, t, are linearly related to the linear and 
angular velocities, v and uj, by a resistance matrix 



K C 
C ti 







(3) 



Here r] is the liquid viscosity. For a particle of character- 
istic size R the translation K ~ R, coupling C ~ R 2 and 
rotation f2 ~ R 3 matrices depend only on the particle's 
geometry. Above, we have tacitly assumed that the resis- 
tance matrix is expressed with respect to a unique fixed 
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point on the particle called the reaction center which re- 
quires that C is symmetric. The third term in Eq. ([3]) 
describes the gyroscopic effect. 

We assume that the particles are uncharged so that 
the electric field does not exert a force on them, F = 0. 
In this case Eq. gives a linear relation between the 
propulsion velocity v and u) 



(4) 



This relation expresses the so-called propeller effect. Ro- 
tation of chiral particles caused by the external torques 
is accompanied by translational motion. Since the cou- 
pling matrix C has opposite sign for particles with op- 
posite chirality (while K and Cl remain the same), it is 
clear that particles of opposite chirality subjected to the 
same torque will move in opposite directions. 

To describe the rotation of the particle we use the body 
axes defined by the unit vectors x.i, x.2, x.3, whose orien- 
tation with respect to the laboratory frame axes x, y, z is 
specified by the three Euler angles <p, 8 and ip |l2j . Their 
evolution equation follows from the balance of angular 
momentum ([3]), 



= QK 1 t x 



(5) 



Here 



sin 8 



sin ip cos ip 

sin 8 cos ip — sin 8 sin ip 
- cos 9 sin ip — cos 8 cos ip sin ( 



(0) 



is the matrix connecting the particle ang ular velocity 
with the derivatives of the Euler angles 12f , and 



(7) 



is the rotational resistance matrix £1 = ri(Cl — Ck^ 1 ^), 
augmented to account for the presence of the intrinsic 
angular momentum L e . 

In general the torques and forces in Eq. (3) consist of 
the deterministic torques of Eq. (2) and random torques 
and forces arising from thermal fluctuations. In this case 
the system can be described by the Fokker-Planck equa- 
tion for the particle distribution function / = f(t, o) [13[ 
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(8) 



and the following expression for the ensemble averaged 
chiral velocity 



V c/l = v 



l/u 



dt d 6 ob(kTKf + fHU). 



(9) 



Here T is the temperature, k is the Boltzmami constant, 



U 



d • E M H, b = -ifi- 1 C(if - Cn^C)' 1 , o 



denotes the particle orientation (specified by the Euler 
angles), and 1Z is a vector derivative operator ToL 14 1. 
whose representation in the body frame is 
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cos ip — cot 8 sin tp 
- sin ip — cot 8 cos ip 

s 'o 1 



sin H 

COS ip 



1 



These equations describe diffusion in the particle orienta- 
tion space and the translational motion associated with 
its rotation relative to the liquid (we assume that the 
spatial distribution of particles is uniform). 

The strength of thermal fluctuations is characterized 
by the dimensionless parameters MH/kT, and dE/kT. 
We consider an electric field of frequency u, E(i) = 
Esin(27wt). In this case another important parameter 
is the ratio of the frequency v to the rotational equilibra- 
tion rate, isrjR 3 /kT. 

For MH/kT, dE/kT <C 1 Eqs. (09]) can be solved by 
perturbation theory. At low frequencies, vr]R 3 /kT <C 1, 
and to second order in dE and first order in MH we 
obtain 



02 



x Rv 



kT 



dE 
kf 



MH 



(10) 



where \ ~ K~ 1 C/R is a dimensionless measure of the 
particle chirality. Note that 03 arises only at fourth order 
in the perturbation theory, leading to the estimate 03 ~ 
XisRvvR 3 (dEMH) 2 1 (kT) 5 . 

In the opposite regime, MH/ kT, dE/kT, vr/R 3 / kT > 
1 , thermal fluctuations may be neglected and the motion 
of particles is described by the deterministic equations 
@, and ((J). In this case, under the influence of 
E and H the particle orientation changes periodically 
with time and traces a cycle in the space of orientations. 
Relation ([JJ shows that propulsion is possible only if this 
cycle is non-self-retracing. This is the analogue of the 
"clam shell" theorem 15( for forced propulsion at low 
Reynolds numbers. It is worth noting that although the 
cycle traced by the linearly polarized electric field is self 
retracing, the cycle traced by body orientations is not. 

At low frequencies where v <C -^fr an< 4 v <C 
(MH) 2 /(dEr/R 3 ), the frequency dependence of the 
propulsion velocity is linear 



0-1,2 ~ X vR{L e /r)R 3 ), 



(11) 



in contrast to the quadratic dependence in the strong 
fluctuation regime (|10|) . The 03 component of veloc- 
ity has a superlinear dependence on the frequency, see 
Eq. (fl2|) below. Therefore at low frequencies the propul- 
sion velocity is confined to the E-H plane. 

The linear frequency dependence of the propulsion ve- 
locity is a feature which often arises in the adiabatic 
regime, where the orientation of the particle corresponds 
to instantaneous equilibrium. This occurs, for example, 
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in the case where the electric field E is circularly polar- 
ized 

For a linearly polarized electric field and constant mag- 
netic field, the origin of the linear dependence on v in 
Eq. (|TT|) is more subtle. In this case the adiabatic ap- 
proximation holds during most of the oscillation period. 
During these intervals the particle orientation is deter- 
mined by the instantaneous values of E(i) and H, and 
the trajectory of the particle is self-retracing. However, 
the adiabatic approximation is violated at time intervals 
where E(i) is approximately zero. During these times the 
equilibrium orientation of the particle is not unique: the 
particle can rotate freely about the axis pointing along 
M, which is aligned with H. Most of the particle's rota- 
tion and propulsion occurs near these instances. As the 
electric field changes sign the particle rotates about M 
by an angle it. The direction of rotation is determined by 
corrections to the adiabatic approximation, which break 
the symmetry between clockwise and anticlockwise ro- 
tations, and lead to a non-self retracing cycle and non- 
vanishing propulsion velocity linear in v. However, if 
L e = 0, upon reversal M — > — M, the direction of the 
rotation and propulsion velocity also reverses. Therefore 
V c /i averaged over realizations of M vanishes. Account- 
ing for the finite value^of L e results in the difference of 
the resistance tensors VL e in states with ±M, and incom- 
plete cancellation of contributions to "V c h from particles 
with different values of M. This leads to Eq. (fTTI) . 

Numerical solutions of the equations of motion, Eq. ([3]) 
confirm this picture. We choose £3 || M, so that the Eulcr 
angle ip(t) corresponds to rotation of the particle about 
M. Typical results of numerical solutions of Eq. Q arc 
presented in Figures [JJ and [2] The time evolution of the 
angle ip exhibits a staircase structure, increasing (on av- 
erage) linearly with t (see the upper panels in FigfT]) . At 
low frequencies there are two steps per oscillation period 
and each step corresponds to a rotation of the particle 
about M by an angle tt, as in the upper left of FigfT] 
In this regime the <j\ and 02 components of the propul- 
sion velocity are linear with respect to the frequency u, 
whereas the 03 component scales as v 1+a with a re 0.25, 
as shown in the lower panel in FigfTJ Propulsion oc- 
curs near the instances t n = n/v, where the electric field 
changes sign. In these intervals the electric field changes 
linearly with time, E(t) re 2nvE(t — t n ). Therefore de- 
viations from adiabaticity depend only on the product 
vE. This gives the following estimate for 03 in the low 
frequency regime, 



cr 3 ~ x v R 



vi]R 3 dE 



(Amy 



(12) 



As the frequency v increases, the character of the motion 
undergoes a series of bifurcations: the step- like character 
of evolution of ip(t) is preserved but the steps become 
separated by several oscillation periods (see the upper 
right panel in FigUJ. This leads to cusp- like features in 
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FIG. 1: (color online) Upper left: time dependence of the 
electric field and the Euler angles at low frequencies. The an- 
gle tjj increases in a stepwise fashion. Upper right: breakdown 
of the staircase at higher frequencies. Lower: log log plot of 
the averaged-over-M chiral velocity versus frequency. The fre- 
quency exponent of o\ and 02 (overlapping straight lines) is 
1 at low frequencies and that of 03 is (approximately) 5/4 as 
discussed above Eq. <|12p . For higher frequencies the motion 
enters a non-adiabatic regime which changes the values of the 
aforementioned exponents, a fact reflected by cusp-like fea- 
tures in the frequency dependence of the propulsion velocity. 



the frequency dependence of the propulsion velocity (see 
the lower panel in Fig. [JJ). 

When L e = 0, particles with opposite values of M ro- 
tate and move in opposite directions in the H-E plane, 
leading to a vanishing planar displacement upon averag- 
ing over the directions of M. In contrast, when L e ^ 
the in-plane displacement averaged over M is finite. This 
is illustrated in Fig.f2]for the case when H and E are per- 
pendicular to each other. 

We now estimate the magnitude of the above effect. 
Assuming that a single domain ferromagnetic particle is 
roughly spherical, we get an estimate for the dimension- 
less parameter characterizing the magnitude of the gyro- 
scope effect, 



Le 



shn 
6rj ' 



(13) 



independent of the particle size. Here n is the volume 
density of magnetic atoms, h is Planck's constant, and 
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FIG. 2: (color online) Particle trajectories in the laboratory 
frame for H and E directed along the x and y axes respec- 
tively. The motion of particles with opposite values of M are 
shown in blue and green. The trajectories originate from the 
same point and move on average in the positive (blue) and 
negative (green) x-direction. The trajectories in the back- 
ground are calculated with L e = 0. In this case the displace- 
ments for opposite M are opposite. The trajectories in the 
foreground are computed with L e 7^ 0. In this case the magni- 
tudes of displacements for opposite values of M are different. 



s is the spin per atom (in units of H), which in different 
materials can lie in the range 1 — 10~ 2 . Using the viscosity 
of water at normal conditions, 77 ~ 10 _2 g/cm-s, and n ~ 
10 23 cm -3 we get L e /Q ~ (s/6) x 1(T 2 . The estimate 
Eq. ([TT]) holds provided the inequalities v < dEj (?/P 3 ), 
v < (MH) 2 /{i-]R 3 dE) and dE, MH > kT are satisfied. 
In modern experiments electric fields in excess of 10 6 
V/m at frequencies 1 MHz have been realized in aqueous 
solutions. Estimating d ~ eR we see that the required 
inequalities are satisfied for a particle size R ~ 100 nm. 
The magnetic restriction, MH/kT ^> 1, is satisfied even 
in weak magnetic fields for a ferromagnetic particle of 
this size. Assuming that the degree of chirality is x ~ 0.1, 
H ~ 10 Gauss and d ~ 10 4 P and using Eqs. (JTTJ) and fTS]) 
with the aforementioned electric fields and frequencies we 
get the estimates 

(71,2 ~ (0.1 — 10) /xm/s, (T 3 ~lmm/s, (14) 

which show that the effect is detectable. 

Finally, we point out the existence of another class of 



effects that are generically related to the one discussed 
above. These are realized by replacing the linearly po- 
larized a.c. electric field E with either a gradient of tem- 
perature VT, pressure VP, or an oscillating magnetic 
field H = Hh, where H = a sin(27rz4) . To describe these 
effects on the phenomenological level one should make 
the following changes in Eq. (JT|) : e -> VT, e -> VP, 
and e — >• h respectively. Calculating the corresponding 
coefficients which are analogues of the <?i,2,3 in Eq. (JTJ) is 
beyond the scope of the present article. 
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